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Abstract 

Many systems where interactions compete with each other or with constraints 
are well described by a model first introduced by Brazovskii. Such systems in- 
clude block copolymers, alloys with modulated phases, Rayleigh-Benard Cells 
and type-I superconductors. The hallmark of this model is that the fluctu- 
ation spectrum is isotropic and has a minimum at a nonzero wave vector 
represented by the surface of a d-dimensional hyper-sphere. It was shown 
by Brazovskii that the fluctuations change the free energy structure from a 
4> to a (p e form with the disordered state metastable for all quench depths. 
The transition from the disordered to the periodic, lamellar structure changes 
from second order to first order and suggests that the dynamics is governed 
by nucleation. Using numerical simulations we have confirmed that the equi- 
librium free energy function is indeed of a 4> e form. A study of the dynamics, 
however, shows that, following a deep quench, the dynamics is described by 
unstable growth rather than nucleation. A dynamical calculation, based on a 
generalization of the Brazovskii calculations shows that the disordered state 
can remain unstable for a long time following the quench. 
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I. INTRODUCTION 



Kinetics of growth can be broadly classified into two categories; nucleation and spin- 
odal decomposition (or continuous ordering). The former applies to situations where the 
initial state is metastable and the latter to where the initial state is unstable. The distinc- 
tion becomes unclear near the limit of metastability and in systems where the concept of 
metastability itself is ill-defined. A class of systems where the definition of metastability 
becomes ambiguous is one where there is a fluctuation-driven first-order phase transition. A 
notable example of a system undergoing such a transition is a symmetric di-block copolymer. 
These systems undergo microphase separation and there is a temperature-driven transition 
from a uniform phase to a lamellar phase Another system which is governed by the 
same dynamical equation as block copolymers, and so exhibits a similar transition, is a 
Rayleigh-Benard Cell §. 

A theoretical model describing these transitions was proposed by Brazovskii in 1975 [[3]. 
Di-block copolymers were shown to belong to the Brazovskii class by Leibler 0] and the 
theoretical predictions regarding the equilibrium nature of the ordering transition have been 
experimentally verified |J. It was shown by Brazovskii, within a self-consistent Hartree 
approximation, that the fluctuations destroy the mean-field instability and lead to a first- 
order phase transition Theories of nucleation and growth have been constructed based 
on the idea that the static, Brazovskii-renormalized theory can provide an effective potential 
for a stochastic Langevin equation HQ]. In this paper we examine the validity of this 
description by using numerical simulations to study the relaxational dynamics of a model 
described by the Brazovskii Hamiltonian and comparing the results to the predictions of the 
"static" nucleation theories and to the predictions of a dynamically renormalized theory || . 
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II. THE EQUATION OF MOTION AND NUMERICAL SIMULATIONS 



Our starting point is the same as in Swift and Hohenberg J7| and describes a system with 
model A dynamics || and a Brazovskii Hamiltonian such as would apply to the diblock 
copolymers ||. The Brazovskii Hamiltonian is characterized by a fluctuation spectrum 
whose maximum occurs at a non-zero wave vector |q| = q and can be represented by a 
hypersphere in d— dimensions. The full form of the Hamiltonian is 

H = J rfx[i0(V 2 + qlf<P + \r<? + 1A0 4 ] (1) 



The dynamics is taken to be relaxational and so the equation of motion is given by the 
Langevin Equation. 



where M is the mobility (which sets the time scale for the problem) and rj is Gaussian noise 
(< v >= 0, < T](x, t) 7](x', t') >= 6(x - x')5{t -t') ). 

The stochastic Langevin equation derived from the Hamiltonian differs from the usual 
Ginzburg-Landau description |J because of the appearance of the unusual gradient term. 
This dynamical equation is usually referred to as the Swift-Hohenberg equation and falls 



into the Type I s classification of Cross and Hohenberg |Tl| . In this classification the system 
is unstable to a static, spatially periodic structure. 
The complete equation of motion is; 

^ = -M(g o 2 V 2 + VV + (4 + r)4> + + r, (3) 



In Eq. (0), the coupling constant, A has been rescaled by the noise strength. For systems 
where the noise strength is small, such as Rayleigh-Benard convection [110], the effective 
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coupling constant is also small. In di-block copolymers, the noise strength is of the order of 
k B T. 

For numerical calculations, the Langevin equation must be approximated by a discrete 
equation. 

0(t + At) = 0(t) - MAt{£<f) + £ 2 <p + (r + q 4 )<f) + ^0 3 ) + (MAt/(Ax) 3 ) 1/2 V (4) 

£ here is the discrete Laplacian. This can take on several forms though for this work we use 
the simplest form: 44 rs £<p = ~^~s( x + Ax) + <p(x — Ax) — 2(f)(x)) where the sum is over 
the dimensions of the lattice. Other choices which include next nearest or more complicated 



neighborhoods are possible ||12|| . This is important if isotropy is of concern however the 
effect is small and can be ignored. The scaling of the noise in Eq. [| reflects the fact that as 
the cell size or time steps become smaller (larger), the possible fluctuations become larger 
(smaller). 

For these simulations we set qo = 1 and choose the lattice spacing such that one lamellar 
spacing spans six lattice sites (Ax = 6/27r) For most of this work the overall lattice size is 
60 3 . The mobility, M, is set to one. In order for the simulations to be numerically stable 
the time scale must be smaller then some stability time, At ~ 1/Ax 4 . To satisfy this the 
timescale is chosen to be At = 0.01 and measurements are taken at intervals of 100 * At or 
longer. 

Previous theories of nucleation have analyzed the above equation under the approxima- 
tion that the effects of fluctuations can be incorporated into a static renormalized free energy 
Fr which replaces the bare free energy F in Eq. fl2|) ||[7|]- We will compare our numerical 
results to the Brazovskii predictions for the static renormalized parameters and show that 
the static scenario is beautifully borne out by the simulations in three dimensions. We 
will then analyze the time-dependence of the structure factor as observed in the numerical 
simulations and show that this time dependence is inconsistent with the dynamical picture 
based on a static renormalized free energy. A theory based on a dynamical renormalization 
of Eq. (§) can provide a qualitatively correct description of the numerical simulations. 



III. STATICS 



A. Brazovskii Theory 

Brazovski's treatment of the model within the Hartree approximation can be restated in 
terms of an expansion of the thermodynamic potential r(</>), the generating functional for 
the vertex functions [0,[L4[]. This approach has been described in detail by Fredrickson and 



Binder M, and within the Hartree approximation leads to a renormalized mass term (r) in 
Eq.([TJ). The diagrams up to one loop are shown in Fig. [IJ, and the mass renormalization 
relation is given by; 



r + {q 2 - ql? = r + (g 2 - g 2 ) 2 + £ -^(r + ( 9 2 - g 2 ) 2 )" 1 (5) 



2 2\ 2 _ ,2 2*2 i _\ f (- i 2 

2 J (2n) d 

This approximation is made self consistent by replacing the bare parameter in the inte- 
grand by the renormalized parameter. Essentially the bare propagator on the loop in Fig. 
(0) is replaced by renormalized propagator. The leads to the Hartree result. 

A f dq 



r = T + 2 J (ip (r + {q2 - ql)2yl = T + aX/V¥ (6) 



where a is proportional to the surface area of a d dimensional sphere of unit radius. Ac- 
cording to Brazovski [^J this approximation is good only for A~ 6 << 1. 

The interesting point about Eq. (|j) is that even for negative values of the bare parameter, 
r, the renormalized parameter, r, is positive. This implies that the disordered phase is 
always metastable. Brazovskii went on to show that the bare coupling parameter A gets 
renormalized to a negative value leading to a 6 theory and the possibility of a first-order 
phase transition 0. 



B. Computational Results. 

The numerical simulations can measure the static structure factor in the disordered phase 
which is predicted by the Brazovskii theory j3]|| to be, 
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S{q)- 1 = r + - qlf (7) 

where r is the renormalized control parameter 0,0]. Thus, S^o) 1 is just r. The renor- 
malized mass can, therefore, be measured by monitoring the peak of the structure factor. 
The Hartree calculations predict r to be positive and asymptotically approaching zero as 
t — ► — oo. Experiments in symmetric di-block copolymers have verified that the behavior of 
S(qo) is consistent with the Brazovskii predictions and is very different from the mean-field 
prediction ( S(q ) = r _1 ). 

Fig. |2| shows results for r obtained from our simulations for different values of the coupling 
constant, A. The Gaussian case (A = 0) is also shown for comparison. The system was run for 
1000 time steps (each time step being 100* At) and 100 samples were taken. The data points 
plotted are the average of the samples while the error bars represent the standard deviation. 
For large positive values of the bare control parameter, r, the fourth order term becomes 
less important and r approaches the Gaussian value. This is seen clearly for small values of 
A. As r approaches zero, the measured values of r deviate from the Gaussian case and stay 
positive for all r. Some care was taken to normalize the values presented in Fig. 0. For the 
Gaussian case, S(qo) _1 should be linearly related to r and the slope of the line should be one. 
In calculating S(qo) several normalizations are needed including the normalization due to 
the Fourier Transform (FT) and the normalization due to circular averaging. While the FT 
normalization is just related to the system size, the normalization due to circular averaging 
is more complicated to calculate. Instead of calculating these normalizations directly, the 
slope of the raw result from the Gaussian case was used to normalize all of the data. Another 
concern is that for an infinite continuous system this line should have an x-intercept at zero; 
however, for the finite systems on a grid used in the simulations the x-intercept is slightly 
negative. To account for this, all of the data is shifted over by an amount equal to that 
intercept. The importance of this shift will become apparent in the next paragraph when a 
scaling is applied and the values of the results around zero are magnified. A small negative 
(positive) value becomes a much larger negative (positive) value and so shifts from negative 
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to positive will become important. 

In reference the authors show that, within the Hartree approximation, the r(r) curves 
for different X's could be described by a single functional form in terms of scaled variables 
t* and r*; 

r* = r * (aA)" 2/3 (8) 
T * = T * (a\)- 2/3 

Figure ^| shows plots of r* against r*obtained from simulations for the three different val- 
ues of A. The curves are seen to scale quite well but the scaled curve falls significantly below 
the theoretical prediction |3|||0] shown as the thick black line. The theoretical prediction 
that the value of r*(and r) never becomes negative and is asymptotic to zero is, however, 
clearly borne out by the simulations. The deviation from the theoretical line could be a 
system size effect, however, results for larger and smaller systems are consistent with the 
data shown in Fig.|3]. For the Gaussian case r* = r* and for large positive values of the bare 
parameter both the theory and simulations approach this. Since the theory is larger then the 
simulations and both are above the Gaussian, our simulations suggest that the contributions 
to the two point function from diagrams not included in the Hartree approximation serve 
to lower the overall correction to mean field theory. Another property that the simulation 
data exhibits in Fig. ^ which is not predicted by the theory is that the curves for different 
values of A diverge from each other at negative values of r, i.e. the scaling is not perfect. 
Smaller values of A lie closer to the theory as expected. The computations were done using 
values of A ~ and larger while the theory is valid only for A ~ 10~ 6 so the theory's 
scaling predictions seem to be far more robust then expected. 

To further explore the nature of the phase transition in this model, we can compare scaled 
values of S(qo) obtained from a hot (random) and cold (purely modulated order) initial 
configurations. Fig. |] shows just such a plot. As in the previous plot, the theory is shown 
as a thick black line. For each value of A there are two lines presented: one for a disordered 
start and one for an ordered start. The averaging was done as described above. For positive 
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and small negative values of r the results from the disordered start and the ordered start 
are nearly the same. For r* below some r* the two differ with the ordered phase having a 
higher peak. It should be pointed out that for these values of r* the disordered state is not 
in equilibrium (metastable or stable) and there is a slow but definite time evolution of the 
structure factor over the time period in which the averages are taken. This time-evolution 
will be discussed in more detail in the section on dynamics. For now the disordered start 
points are included in the plot simply as a comparison to see where the ordered phase melts. 
The value where the two data sets diverge can be considered an estimate for the limit of 
stability (spinodal) of the lamellar phase. As a consistency check the average value of the 
wave vector was measured. For values of r above t s the wave vector is small and points in a 
random direction, while for lower values of r the average wave vector is large and and points 
along the direction that the system was prepared in. For all three sets of data the lamellar 
spinodal lies in the range of —2.1 < r* < —1.9. This is consistent with the prediction from 
Hartree theory that the first-order transition to the lamellar phase occurs at a lower value 
of r*, r t * = —2.74 J7J. That is to say that our measured value of the lamellar spinodal 
is above the predicted transition temperature as it should be. For values of r* below r*, 
however, systems prepared in the disordered phase always evolve towards the ordered phase. 
It is possible that the lamellar ordered phase is enhanced by the boundary conditions - the 
lattice size is always chosen to be commensurate with the lamellar wavelength - however, it 
should be noted that lamella form in many possible directions where the system size is not 
commensurate with the lamellar wavelength. Discussion in the dynamics sections should 
shed some light on this issue. 



IV. DYNAMICS 



S 



A. Relaxational Dynamics Based on Static Renormalized Parameters 



Previous analysis of nucleation and metastability, in the context of fluctuation-driven 
phase transitions, have been based on a Langevin equation with the force obtained from 
the static, Hartree-renormalized free energy function |],|7||. Fredrickson and Binder || used 
this approach to compute the nucleation barrier and the completion rate of nucleation and 
growth in di-block copolymers. The interfacial tension was found to be small leading to 
a small nucleation barrier and rapid nucleation for deep quenches. In this picture, there 
is no essential difference between the kinetics of a fluctuation-driven first-order transition 
and a weak-first order transition. The results of our simulations suggest a very different 
scenario for the growth of the lamellar structures in a Brazovskii model. The shape of the 
nucleating droplet was more carefully analyzed by Swift and Hohenberg by taking into 
account spatial inhomogeneities in the effective free energy function. Their analysis relied 
on constructing a coarse-grained free energy functional. The coarse graining was based on a 
momentum-shell renormalization idea where fluctuations with momenta far away from the 
shell defined by |q| = qo are successively integrated out. This analysis led to non-spherical 
droplets and were consistent with the picture of spinodal nucleation that had been obtained 
earlier ]15| . 



As pointed out in the work of Swift and Hohenberg [|7| , a complete theory of the dynamics 
of fluctuation-driven first order transitions would have to be based on a coarse-graining of 
the full dynamics as expressed by the original Langevin equation with the bare Brazovskii 
Hamiltonian. In this paper, we compare the results of our numerical simulations to the 
predictions of a Hartree-renormalization of the full dynamical equation as expressed in Eq. 
(§). Our emphasis has been on understanding the nature of the dynamics of the metastable 
phase and we have not analyzed, in any detail, the spatial structures associated with the 
growth process. 
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B. Dynamical Renormalization and Simulation Results 



In order to systematically analyze the effects of fluctuations on the kinetics of growth of 
the lamellar phase, perturbative techniques analogous to the static Hartree approximation 
have to be applied to the Langevin equation. This is most conveniently done through the 
dynamical-action formalism ||14|| . The application of this method to the S-H equation is 
outlined by Ignatiev et. al. ||. In the dynamical-action formalism, the average value of an 
operator 0(<f)) over the noise history is rewritten as a functional integral; 

(O{0(x,t)}) = jD<f>exp(-S[<f>]) , (9) 



where S[(f>] is the dynamical action [14 

As in the static case, the calculation of dynamical correlation functions is most conve- 



niently formulated through the construction of a generating function |14|,|T3J. To establish the 
closest analogy with the static calculations, it is useful to work with the generator of vertex 
functions, r^ n [0], and establish a diagrammatic expansion which is the exact analog of 
the diagrams retained by Brazovskii in the static calculation. The simulation results demon- 
strate that the static properties of the S-H equation closely follow the Brazovskii predictions 
and, therefore, it seems appropriate to apply this approximation scheme to the dynamics. 
The correspondence between statics and dynamics becomes particularly transparent in a 
super-field formulation |I| which shows that the dynamical perturbation theory in terms 
of the super fields has exactly the same structure as the static perturbation theory except 
for the appearance of a very different "kinetic" term which leads to a bare propagator that 
is distinct from the static theory. The super-field correlation functions can be calculated 
by constructing diagrams as in the static theory but replacing the static bare propagator 
by the appropriate dynamical one |]IE|| . The super-field correlation functions encode all the 
dynamical correlations of the field <fi and the latter can be extracted from well-defined re- 
lations [|TIJ. The details of this calculation and the results will be described in a separate 
publication ||. In this paper we present the main results which can be compared to the 
numerical simulations. 
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The two correlation functions which appear in the dynamical study are 



G(t,f) =< 4>{t)r]{t') > 



and 



C(M') =< <t>{t)<t>{t') > 



Both of these can be obtained from a single super-field correlation function (Q(l, 2)) which 
is related to T dyn through 



where $'s are the super fields and the derivative is evaluated with the source term set to 
zero. Applying this formalism to Eq. (|]), leads to an expression for Td yn which involves 
two-point correlation functions of the fluctuations and is a natural extension of the static 
T(0) to correlation functions involving time. The expression, however, does not lend itself 
to easy analysis except in two cases, early times when a variant of linear theory can be 
applied and late times when the system is stationary. The late time dynamics involves the 
analysis of the nonlinear terms and will be discussed in the concluding section. The early 
time dynamics is where one is justified in retaining only quadratic terms in the renormalized 
Tdyn- In this limit we obtain the following equation for the equal-time correlation function 
C q (t,t) =< (p q (t)(p- q (t) > which is the structure factor that is monitored in the simulations: 



The mass parameter, r(t), which is now time-dependent, is renormalized in the dynamic 
theory using the same approximation as in the static theory. The diagrams used are given 



Q (1,2) 




(10) 




(11) 



in fig. [7]. These again are only to one loop. The mass term then becomes 



r(t) + {q 2 - qlf 



r+(q 2 ~ qlf + 
A r dq 1 



(12) 
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Again, the approximation is made self consistent by replacing the r with r in the integrand. 

r(t) -r - A f e Mz3EM ( 13 ) 
T[t) ~ Teq 2 J J^f D(q) [ } 

D(q) here is the renormalized propagator, D(q) = r eq + (q 2 — q%) 2 where r eq denotes the 
static renormalized mass parameter which is the solution to eq. [?]. 

In Eq. [13], at long times the integrand in the second term becomes small and r(t) ap- 
proaches r eq . which is confirmed in the section on statics above. At time zero the subtraction 
of the second term is just a subtraction of the Hartree level correction factor to the bare 
parameter introduced earlier Eq. ((^ and, therefore, r(t = 0) is just the bare parameter r, 
while as t —>■ oo, r(i) approaches the renormalized, equilibrium value. 

This scenario is confirmed by the simulations. The simulation results discussed in Section 
III B confirm that the equilibrium value, r eq , is consistent with the Hartree prediction. 
Analyzing the early stage dynamics can verify the value of r(t) at short times. Figure 
|5] shows the growth of S (go) as a function of time for various values of the bare control 
parameter, r. To average over noise, five independent runs where taken for each value of 
the control parameter. If r(t) is time-independent, the early-stage evolution (linear theory) 
is described by W7 



S(q, t) = C q {t, t) = S(q, 0) eM-D(q)t} + -^r (1 - exp{-D{q)t}) (14) 

where D(q) is D(q) = r + (g 2 — q^) 2 and r is the time-independent value of r. At the peak 
of the structure factor, q = q , and D(q Q ) is just r , the value of which can be estimated by 
fitting the simulation results shown in Fig. |5] to [14|. The results of such fits are presented 
in Fig. H This figure shows that r$ and r are linearly related with a slope which is nearly 
one. This implies that for short times, the dynamically renormalized parameter is linearly 
related to the bare parameter and is not renormalized to positive values for negative values 
of the bare parameter. The very early-stage dynamics, therefore, is characteristic of a system 
exhibiting unstable growth. 

These fits must be considered with some care. In the case where r is not time dependent, 
the linear theory describes the system only for short times, that is times on the order of 
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the natural time scale, r 1 [See for example fT8| . For the fastest growing mode, when the 



data is fit to different time scales the results obtained are consistent for all time scales up 
to r~ l . For the S-H equation, however, the results obtained are dependent on the time 
range. In particular, as the range of time increases the values obtained for r become less 
negative. This is consistent with r(t) growing as a function of time and the results of the 
fits merely represent some average value of r(t) for the time range involved. If the system 
has been quenched to a negative value of r, Eq. [H| indicates that the early-time evolution 
will exhibit unstable growth. As time evolves the second term in Eq. (|13|) decreases and 
the value of r(t) approaches r eq which is always positive. Though the integral in Eq. [L3| is 
hard to evaluate analytically, numerics can provide some insight into its behavior. Figure 
|8] shows a numerical evaluation of the time evolution of r(t) for various values of r eq and 
A = 0.01. As confirmed in the statics section above, r(t) will eventually become positive no 
matter how deep the quench is, though the time for this to happen may become very long. 
As long as r(t) is negative the system is expected to undergo unstable growth though, since 
the mass parameter is time-dependent, the time evolution will be more complicated than 
the usual spinodal decomposition scenario. We can define a crossover time, t cross at which 
r(t) changes sign, and an average growth time, t ave . The average growth time is just the 
inverse of r ave defined as 



ter. For shallow quenches t cross is small compared to t ave and the system reaches a well- 
defined metastable equilibrium state which is disordered. If t cross becomes larger then t ave , 
metastability becomes difficult to define. The system takes longer to equilibrate in the 
metastable disordered phase than it takes to grow lamellar structures. We can use the 
condition t cross = t ave to define a crossover temperature r*[ yn . Above this temperature the 
disordered phase quickly becomes locally stable and a nucleation event is needed for the 
formation of the lamellar structures. Below , the system is expected to evolve contin- 




(15) 



In Fig. [], we show plots of t, 



cross 



and t nvP function of the bare control parame- 
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uously towards the lamellar phase without any evidence of metastability of the disordered 
phase. Fig. |9| is in sharp contrast to standard Ginzburg-Landau (G-L) theory. For the time 
dependent G-L equation, the presence of noise suppresses the critical point. In the region 
where the bare parameter is negative but the renormalized parameter is positive, t cross is 
always smaller then t ave and so there is no unstable growth until the true critical point is 
reached. 

For the system under consideration here, different scenarios are possible depending on 
the relationship of this dynamical crossover temperature to the first-order transition tem- 
perature, r* rans , at which the free energy of the lamellar phase becomes lower that that of 
the disordered phase. If r t * ans is higher than , there will be a regime of temperatures 
over which the the system will undergo nucleation and growth. On the other hand if T t * ans 
is lower than r^ yn , then there is no nucleating regime and one would observe only unstable 
growth, albeit of an unusual nature since the free-energy surface is evolving with time. From 
figure [| we can see that t cross is small only for very shallow quenches below the mean-field 
transition. The results of our simulations, presented in the next section, are in qualitative 
agreement with this scenario. We would like to emphasize that the growth dynamics de- 
scribed above are qualitatively different from rapid nucleation in which the system quickly 
reaches the stable equilibrium state. A large value of t cross , on the other hand, indicates a 
type of ergodicity breaking as the system takes a very long time to reach equilibrium. 

All of our numerical results indicate that T* rans predicted from Hartree is lower than 
T% yn and they are remarkably close to each other. We have been unable to come up with a 
simple relationship between these two temperatures and, therefore, can only interpret the 
similarity of the two as a remarkable coincidence. The dynamical crossover is deduced from 
time scales which characterize the evolution of the free-energy surface while the transition 
temperature is deduced from a comparison of the depth of the two wells. It is not clear 
why the two temperatures should be similar in magnitude. It can be argued that there 
is only a single parameter, A, controlling the scale of fluctuations and, therefore, the two 
temperatures should be related, however, there is no obvious argument to suggest that they 
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should be identical. 

The picture emerging from the dynamical renormalization is a natural extension of the 
effect of fluctuations on the static results of the Brazovskii model. Following a quench from a 
relatively high temperature to a temperature where r is negative, the system is in a locally 
unstable region (top of a hill). As the fluctuations grow with time, the non-linear terms 
characterized by A become important and they renormalize the curvature of the hill. This 
scenario is quite different from the usual picture of evolution in an adiabatic potential. 



C. Simulation Results for Late Times 

Further evidence for the dynamical scenario presented in the previous section comes from 
examining the long time evolution of the peak of the structure factor obtained from our 
numerical simulations and correlating that to snap shots of the system as it evolves. Figure 



ITJ is a plot of the amplitude of the structure factor peak as it evolves. For shallow quenches, 
t > —0.04, the peak grows to an equilibrium value and does not evolve further. These values 
are consistent with the equilibrium values reported above in the section on statics. For 
deeper quenches the peak amplitude grows quickly and then appears to saturate; however 
careful examination shows that the value continues to grow very slowly. This is consistent 
with late stage domain growth which as been studied extensively by Elder et. al. for 2-D 



systems |22] . The last feature in the system is a final rise to an equilibrium value. This rise 
is a finite size affect is caused by the majority domain finally taking over the entire system. 
This being the case, the peak would not be expected to grow much since the difference 
between the value before and after this final evolution should only represent the surface area 
between the different domains. 

For a quench depth of r = —0.07 this time evolution can be compared to a series of system 



snapshots shown in figures [TTJ- [15|. These figures represent = isosurfaces which would be 
the boundaries between the different microphases of the system. For early times the system 
appears to be very disordered while at later times domains of ordered lamellar structures 
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begin to appear. Just before the final convulsive growth, at a time t = 4000, there still 
appear to be a few different ordered domains while just after the final growth, at t — 5000, 
only one domain appears to be left in the system. Thus we still have a consistent picture of 
slow but continuous domain growth in the system which is governed by the evolution of the 
2nd order term as it goes from negative values to positive values. 



In the work presented above we have shown that analysis of the static free energy does not 
always provide an adequate description of the fluctuation driven dynamics. Though figure 
|10| shows features which are reminiscent of a first order phase transition, dynamical analysis 
of the renormalized coefficients suggest that a more complicated evolution is taking place. A 
similar situation occurs for the superconducting transition which is also a fluctuation-driven 
first-order transition ||19|| . 

In the present work we have paid close attention to the time evolution of the mass 
parameter and its effect on the evolution of the system at early times. The dynamical 
renormalization results also become tractable at times large compared to t cross when the 
system has reached equilibrium. As pointed out earlier, one enters this regime fairly quickly 
for shallow quenches. At this stage, the higher order terms in Td yn become important. The 
most interesting aspect of the renormalized theory is that the fourth-order term becomes 
non-local in time |§: 



and the integral of the kernel over all times leads to the Brazovskii result for the renor- 
malized fourth-order term. One can derive an "effective" Langevin equation for the fields 
(f) from the renormalized Tdyn, and the above equation implies that this Langevin equation 
has "memory" effects which appear in the non-linear term. This would lead to unusual 
behavior of the equilibrium correlations functions and might provide a distinctive signature 



V. FUTURE WORK 
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of fluctuation-driven first-order transitions. We are in the process of investigating the ef- 
fective equation for numerically. One obvious consequence of the non-local term is that 
the barrier to nucleation is dynamic and it is not valid to think of nucleation as tunneling 
through a fixed barrier. The situation is closer to the problem of quantum tunneling for 
many degrees of freedom where one also finds an effective equation which is non-local in 



time 20 



Another area of interest is the role of defects in the phase transition of the system. 
These will be important for very shallow quenches to just below the transition point and 
we have noted there presence when looking at snapshots of the system. These defects are 
reminiscent of the perforated lamellar phase described in reference [[I]. One way to study 
this is to use a directional order parameter measure introduced by Bray et al. [^TJ to study 
the 2D problem. This order parameter is similar to the order parameter for complex fluids 
and provides directional as well as density information. What it may show is that the 
local gradient terms are enhanced for negative values r suggesting that perhaps there is a 
disordered to nematic transition in this region which others have suggested for 2D [2~2 . 

We would like to close with a discussion of possible experimental systems where this 
dynamical scenario can be studied experimentally. Since the effective A is extremely small in 
Rayleigh-Benard systems the regime where the disordered state is metastable will probably 
be unaccessible. In di-block copolymers, however, there should be a significant temperature 
range, below the mean-field transition, where the disordered state is locally stable and a 
study of the equilibrium, two-time correlation function should reveal the existence of the 
non-linear memory term. It might also be possible to observe unusual nucleation if the 
system parameters allow for r t * ans to be higher than r*[ yn . 
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FIGURES 
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FIG. 1. Diagrams used for static renormalization of the mass parameter in the free energy. The 
Hartree approximation only uses diagrams up to order A and then replaces the bare parameter in 
the integrand with the renormalized one to solve the equations self consistently. 
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FIG. 2. The renormalized control parameter, r, against the bare control parameter, r. The 
different data sets are for different values of A. The Gaussian case, A = 0, which was used to set 
the scale, is also plotted for comparison. The inset shows the region near the mean field transition. 
Note that the renormalized parameter does not cross zero here. 
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FIG. 3. Scaling of r(r) using the form deduced from Hartree theory (c/ text). The data 
collapses well except at the deepest quenches. The scaling curve predicted by the Hartree approx- 
imation is shown for comparison. It is in fairly good agreement with the data considering that the 
values of A used in the simulations are four order of magnitude larger then that for which the the 
approximation is expected to be valid. 
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FIG. 4. Temperature dependence of the average intensity at the peak of the structure factor. 
To show the transition, results from two different starting configurations are compared. The shapes 
of the symbols denote different values of A, whereas the shading distinguishes between the starting 
configurations. The closed symbols represent systems that were started in a disordered state, 
simulating a hot start, while the open symbols are for systems which were started with an ordered 
lamellar structure already present. At high temperatures the two starting conditions give nearly 
the same results while at low temperatures the systems prepared in the ordered state have higher 
peak intensities implying that the systems prepared in the disordered states are in a metastable 
state or not in equilibrium (c/text). 
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FIG. 5. The time evolution of the structure factor peak for the earliest times. Different sets 
are for different quench depths as indicated on the graph. All quenches to below the mean field 
transition temperature initially show unstable growth. 
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FIG. 6. Growth time as calculated from fits to a linear theory plotted against the quench depth 
below the mean field transition. For short times the system should be described by linear theory 
which predicts unstable growth for all quenches below the mean field transition. 
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FIG. 7. As with the static case, diagrams only up to order A are used and then the bare 
parameter is replaced by the renormalized parameter to give a self consistent solution. 




FIG. 8. The time evolution of the mass parameter as calculated by dynamic renormalization. 
r(t) eventually becomes positive for all quenches below the mean field transition but the time for 
that cross over increases dramatically for larger quenches. At large times r(t) approaches r eq . 
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FIG. 9. The calculated crossover time plotted against quench depth. Also plotted is a char- 
acteristic growth time for the lamellar structure. These results have been scaled using equation 
H so that they are A- independent. When this time becomes much less then the crossover time 
then the lamella has enough time to grow. Notice that this occurs just above or at the transition 
temperature, r 4mns . 
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FIG. 10. The evolution of the structure factor peak as a function of time for different quench 
depths. For shallow quenches above the real transition temperature S(qo) grows quickly to its 
equilibrium value. For quenches to below the transition temperature the system grows quickly at 
early times and then reaches a late stage regime where it evolves by domain growth and S(qo) 
grows as a small power law in time. The highlighted line is a quench to r = —0.07. Snapshots of 
this system are presented in the next figure. 
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FIG. 11. Time = 1000 




FIG. 12. Time = 2000 



FIG. 13. Time = 3000 
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FIG. 14. Time = 4000 



FIG. 15. Time = 5000 
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